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Abstract.  Collisional-radiative  model  for  atomic  hydrogen  has  been  coupled  selfconsistently  with  the  Boltzmann 
equation  for  election  kinetics  to  study  the  role  of  nonequilibrium  electron  distributions  on  the  energy  level  kinetics. 
In  particular,  the  role  of  second  kind  and  electron-electron  collisions  in  affecting  both  time  dependent  electron  energy 
distributions  and  excited  state  population  densities  has  been  investigated. 


INTRODUCTION 


The  analysis  of  the  emission  spectra  of  a  plasma  is  a  powerful  tool  to  investigate  plasma 
composition  and  temperature.  Quantitative  analysis  is  based  on  the  proportionality  law  that  relate  the 
amplitude  of  a  line  and  the  population  density  of  the  emitting  level.  New  techniques,  such  as  LIBS 
(Laser  Induced  Breakdown  Spectroscopy),  have  renewed  the  interest  in  emission  spectroscopy  and 
collisional-radiative  (c-r)  models  for  atomic  systems.  In  spatial  homogeneous  plasmas,  the  population 
of  excited  levels  is  determined  by  collisional  processes  (mainly  electron-atom,  e-A)  and  radiative  decay. 
Collisional-radiative  models  calculate  level  distributions  solving  the  system  of  master  equations  (one 
equation  for  any  level)  where  the  kinetic  terms  are  due  to  e-A  collisions,  radiative  decay  and  absorption. 

The  classical  approach  to  c-r  models  is  the  QSS  (quasi-stationary  state)  approximation  introduced 
by  Bates  et  al.  [1].  This  model  is  based  on  the  assumption  that  the  relaxation  times  of  excited  levels  are 


order  of  magnitude  shorter  than  ionization/recombination  ones  (-^  =  0  Vi  ^  0,— =  -^£- 

3t  dt  3t 


*0).  In 


these  conditions  the  level  distributions  are  stationary  with  respect  to  electron  and  ion  populations, 
simplifying  the  solution  of  the  problem. 

Soon  after,  however,  it  was  recognized  the  importance  of  a  time  dependent  integration  of  kinetic 
equations  which  are  at  the  basis  of  the  c-r  models  for  situations  characterized  by  very  short 
macroscopic  relaxation  times.  Time  dependent  solution  of  c-r  models  have  been  presented  by 
Cacciatore  et  al.  [2-5]  for  atomic  H,  O  and  N  plasmas  many  years  ago  and  more  recently  these  plasmas 
[6-8]  have  been  again  studied  for  explaining  new  experimental  situations  (e.g.  laser-plasmas  interaction) 
in  which  a  time  dependent  solution  of  c-r  models  seems  to  be  necessary.  In  these  papers  [2-5],  rate 
coefficients  of  e-A  collisions  are  calculated  approximating  the  eedf  with  the  Maxwell  distribution  at  a 
fixed  temperature  (free  electron  kinetics  have  been  neglected).  This  assumption  is  realistic  only  if  the 
ionization  degree  is  high  enough  so  that  electron-electron  Coulomb  collisions  (e-e)  dominates  the  free 
electron  kinetics. 
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To  applicate  c-r  models  to  the  region  in  which  the  actual  eedf  differs  appreciably  from  the 
Maxwellian  form,  in  the  present  work  we  have  coupled  selfconsistently  a  collisional-radiative  model  for 
hydrogen  atoms  with  the  Boltzmann  equation  for  free  electron  kinetics,  calculating  simultaneously  level 
distributions  and  eedf.  This  model  allows  to  investigate  departures  from  equilibrium  (Maxwell  and 
Boltzmann)  distributions  during  the  thermal  relaxation  of  the  hydrogen  plasma. 

The  aim  of  this  paper  is  to  analyze  the  conditions  in  which  the  electron  energy  distributions  and  the 
excited  state  populations  differ  from  Boltzmann  and  Maxwell  ones.  Moreover,  we  are  interestedin 
studying  the  role  of  nonequilibrium  electron  distributions  on  the  level  kinetics.  Finally,  we  intend  to 
understand  the  influence  of  each  elementary  processes  considered  on  both  the  time  dependent  electron 
energy  distribution  and  the  excited  state  population. 


METHOD  OF  CALCULATION 

The  system  studied  in  this  work  is  a  plasma  formed  by  atomic  hydrogen  (H),  protons  (H+)  and 
electrons  (e').  We  assume  the  plasma  is  homogeneous  and  neutral  (ne=n+).  The  levels  are  grouped  by 
principal  quantum  number  up  to  a  maximum  value  of  25. 

The  plasma  has  been  considered  optically  thin  (photon  absorption  and  photoionization  have  been 
neglected).  The  collisional-radiative  model  consists  in  a  set  of  master  equations,  one  for  ion  (n+)  and 
electron  (ne)  molar  fraction  and  one  for  each  H  level  (nj) 
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where  Ay  are  the  Einstein  coefficients  of  the  radiative  decay  (process  3),  Rj  the  rates  of  radiative 
recombination  (process  4),  KT  and  K®.,  with  i<j,  are  respectively  the  rates  of  electron  impact  excitation 

(forward)  and  deexcitation  (backward)  (process  5),  S®  and  are  the  rates  of  electron  impact  ionization 
and  three-body  recombination  (process  6) 
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where  Ey  and  £,  (eq.  7)  represent  the  threshold  energy  of  the  transition  i— >j  and  the  ionization  energy 
from  the  i-th  level.  Atom-atom  and  ion-atom  collisions  have  been  neglected  assuming  that  the  collision 
frequency  of  electron-atom  is  larger  as  both  mean  electron  speed  and  electron  impact  cross  sections  are 
higher  than  heavy  particle  ones. 

Rate  coefficient  for  electron-atom  collisions  are  calculated  from  eedf  using  the  equation 

Ke  =  J“  f  (£)o(£)v(E)d£  (8) 

where  f(£)  represents  the  eedf,  o(£)  the  cross  section  of  the  transition  between  atomic  levels  i  and  j,  v(£) 
the  velocity  of  the  electron  of  kinetic  energy  £  and  Et  the  threshold  energy  of  the  process.  The  inelastic 
cross  sections  are  taken  from  databases  available  on  internet  [9],  as  well  as  the  Einstein  coefficients. 

Cross  sections  for  superelastic  collisions  and  three-body  recombination  are  calculated  using  the 
detailed  balance  principle  [10,11]: 
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where  gj  is  the  statistical  weight  of  the  i-th  level  and  gj+  is  the  statistical  weight  of  the  ion  ground  state, 
h  the  Planck  constant,  me  the  electron  mass,  k  the  Boltzmann  constant  and  Te  the  electron  temperature. 
Radiative  recombination  rate  coefficients  have  been  calculated  according  to  the  work  of  Johnson  [12]. 

To  tintroduce  the  non  equilibrium  free  electron  distributions  we  have  solved  the  Boltzmann  equation 
in  the  homogeneous  and  quasi-isotropic  approximation  developed  by  Rockwood  [11,13]  in  the  time 
dependent  form 


0f(£,t)  _  0Je  djc]  dJe-e  c  c 

- -  — — - — - T - Oin  —  OS  up 

3t  d£  d£  de  (11) 

where  f(£,  t)  is  the  time  dependent  electron  energy  distribution  function  (eedf).  The  J  terms  on  the 
right-hand  side  of  eq.  1 1  are,  respectively,  the  flux  of  electrons  in  the  energy  space  due  to  the  electric 
field  (Jg),  elastic  (Je|)  and  electron-electron  (Je_e)  collisions  and  the  S  terms  consider  the  jump  of 
electrons  in  the  energy  space  due  to  inelastic  (Sjn)  and  superelastic  (Ssup)  collisions.  Ionization  of  atoms 
by  electron  impact  (see  eq.  6)  has  been  considered  in  the  electron  kinetics  as  an  inelastic  process, 
neglecting  the  production  of  secondary  electrons,  while  in  the  master  equations  (eqq.  1,2)  this 
production  has  been  taken  into  account.  The  numerical  solution  of  this  equation  is  widely  described  in 
literature  (see  for  example  [11,  13]). 

The  two  systems  of  equations  (one  for  electron  kinetics  and  the  other  one  for  the  collisional-radiative 
model)  are  strictly  correlated;  rate  coefficients  of  electron  impact  processes  are  calculate  from  the  eedf 
(see  eq.  8)  and  the  Boltzmann  equation  terms  J  and  S  (see  eq.  11)  depend  on  the  level  distribution.  The 
simultaneous  solution  of  electron,  level  and  ionization  kinetics  (self-consistent  model)  is  therefore 
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necessary.  To  calculate  the  time  evolution  of  eedf  and  level  distribution,  we  have  used  the  first  order 
implicit  Euler  method  with  step  adaptive  algorithm  [14]. 


RESULTS  AND  DISCUSSION 

In  this  section  we  will  investigate  the  relaxation  of  the  level  population  and  eedf  in  different 
conditions.  The  initial  state  is  defined  by  fixing  the  gas  temperature  and  pressure  (maintained  constant 
in  all  the  integration  interval),  and  assigning  the  initial  value  of  ionization  degree  (ID),  level  population 
(a  Boltzmann  distribution  with  temperature  Tj)  and  eedf  (a  Maxwell  distribution  with  temperature 
Te).  In  all  the  cases  studied,  the  pressure  and  initial  ionization  degree  have  been  fixed  respectively  to 
the  value  of  P0=l  atm  and  ID0=10'3,  and  the  cases  considered  differ  for  the  value  of  T,  Tj  and  Te.  The 
gas  temperature  does  not  influence  directly  the  level  kinetics  as  atom-atom  and  atom-ion  collisions  have 
been  neglected.  On  the  contrary,  the  relaxation  time  of  the  eedf  due  to  e-H  and  e-H+  elastic  collisions 
depends  on  gas  temperature  and,  as  a  consequence,  indirectly  affects  the  level  distribution  through 
inelastic  and  superelastic  collisions.  We  have  analyzed  three  different  conditions  listed  in  table  1 . 


TABLE  1.  Initial  conditions  for  the  cases  studied 


Label 

T  (K) 

Ti(t=0)(K) 

Te(t=0)(K) 

ID() 

P0(atm) 

1 

1000 

10000 

20000 

lO-3 

1 

2 

1000 

10000 

5000 

10-3 

1 

3 

20000 

10000 

1000 

10-3 

1 

Figure  1  reports  level  population  and  eedf  in  the  case  1  of  table  1  (T<Tj<Te).  The  level  population 
(n;)  has  been  divided  by  the  level  statistical  weight  (gj)  since  for  a  Boltzmann  distribution  the  logarithm 
of  nj/gj  shows  a  linear  dependence  on  the  level  energy. 


FIGURE  1.  Time  evolution  of  level  distribution  (a)  and  eedf  (b)  in  the  case  1  of  table  1. 
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In  this  case  we  have  Te>T  and  therefore,  due  to  electron-atom  and  electron-ion  elastic  collisions,  eedf 
cools  down  to  equilibrate  the  gas  temperature.  It  can  be  observed  that  the  electron  distribution  tail  is 
non-maxwellian  also  at  t<10'9  s,  while  the  electron  temperature  (which  depends  on  the  low  energy  eedf) 
does  not  change.  Because  Te>Tj,  inelastic  collisions  and  electron  impact  ionization  are  the  dominant 
processes  in  eedf  and  level  population  relaxation.  Only  when  the  electron  distribution  is  cold  enough 
(t>10  6  s),  the  effects  of  superelastic  collisions  affect  the  eedf  (see  plateaux  in  fig  1  b). 


-15  -13  -11  -9  -7  -5 

10  10  10  10  10  10 

time  (sec) 


FIGURE  2.  Time  evolution  of  ionization  degree  in  case  1  of  table  1. 

While  the  electron  distribution  is  cooling,  level  distribution  tail  decreases,  deviating  from  the 
Boltzmann  distribution.  For  t<10"9  s,  level  kinetics  is  dominated  by  ionization,  while,  for  t>10"4  s,  the 
spontaneous  emission  for  the  low  excited  levels  prevails  and,  at  the  same  time,  radiative  recombination 
tends  to  repopulate  the  last  excited  levels.  This  behaviour  can  be  confirmed  by  time  evolution  of  the 
ionization  degree  (see  fig  2).  Ionization  degree  increases  in  the  time  interval  10~12-rl0'10  s,  when 
ionization  prevails,  while  it  decreases  for  t>10"5  s  due  to  e-H+  recombination. 

In  case  2  of  table  1  (fig.  3  a,  b),  electron  temperature  is  higher  than  gas  temperature  but  lower  than 
internal  temperature  (T<Te<Ti).  Initial  (t=0)  and  final  (t=°°)  eedfs  are  maxwellian  on  the  contrary 
transient  states  show  nonequilibrium  distributions.  Two  different  trends  are  observed:  for  t<10"9  s,  eedf 
warms  up  while  it  cools  down  for  t>10"9  s.  Firstly  (t<  1 0“ 1 2  s)  two  peaks  appear  (due  to  superelastic 
collisions)  and,  for  longer  time  (t=10  9  s),  elastic  collisions  (e-H,  e-H+,  e-e)  transform  them  in  plateaux, 
resulting  in  the  overpopulation  of  eedf  tails.  For  t>10'9  s,  the  energy  exchange  between  electrons  and 
excited  levels  equilibrates  and  therefore  the  eedf  cools  down  slowly  because  elastic  collisions  push  the 
eedf  towards  a  maxwellian  at  the  gas  temperature.  The  level  distribution  cools  down  as  well  as  does  the 
eedf  and  also  as  a  result  of  the  radiative  decay.  Similarly  to  case  1,  for  sufficiently  long  times,  the  tail  of 
level  distribution  is  overpopulated  with  respect  to  the  Boltzmann  one  when  the  radiative  recombination 
becomes  more  important  than  ionization. 

Very  interesting  results  are  obtained  when  the  initial  electron  temperature  of  the  plasma  is  lower  than 
both  internal  and  gas  temperature  (Te<Tj<T),  as  in  case  3  of  table  1  (fig  4  a,  b).  Again,  initial  and  final 
eedf  are  maxwellian  but  intermediate  distributions  are  very  far  from  equilibrium.  In  particular,  at  short 
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time  (t=10'n  s)  the  distribution  tail  shows  sharp  peaks  because  of  superelastic  collisions  and  at  longer 
times  (t=10~8  s)  plateaux  appear  due  to  synergetic  contribution  between  elastic  and  superelastic 
collisions.  Eedf  tends  to  increase  its  temperature  and,  when  it  is  hot  enough,  the  inelastic  collisions 
prevail,  determing  a  change  in  the  eedf  slope  at  the  threshold  of  the  first  inelastic  process.  This 
behaviour  is  reflected  on  the  level  distribution  (see  fig  4  a),  which  cools  till  t=5  10~8  s,  because  the 
energy  is  transferred  to  electrons  through  superelastic  collisions,  and  later  it  is  populated  by  inelastic 
collisions. 


FIGURE  3.  Time  evolution  of  level  distribution  (a)  and  eedf  (b)  in  the  case  2  of  table  1. 


FIGURE  4.  Time  evolution  of  level  distribution  (a)  and  eedf  (b)  in  the  case  3  of  table  1. 
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To  study  the  role  of  non  equilibrium  electron  distributions  on  level  kinetics,  we  have  performed 
calculations  in  the  same  condition  of  case  3  but  fixing  eedf  to  a  maxwellian  distribution  at  the  initial 
electron  temperature  during  all  the  temporal  evolution.  In  this  approximation,  we  neglect  both  eedf 
deviations  from  equilibrium  and  variations  in  electron  temperature  and  consider  constant  rate  coefficient 
for  electron-atom  collisions  (see  fig  5).  Comparing  fig  4  a  with  fig.  5,  it  can  be  noted  that  level 
population  of  fig  5  undergoes  only  a  cooling  and  level  population  remain  constant  for  10  7  s<t<  1CT8 
s,  while  in  fig  4  a,  in  the  same  time  interval,  level  population  increases  by  3-4  orders  of  magnitude. 
These  results  show  that  time  dependent  selfconsistent  calculations  are  necessary  to  have  an  accurate 
study  of  the  time  evolution  of  level  population. 


level  energy  (eV) 


FIGURE  5.  Time  evolution  of  level  population  in  case  3  fixing  eedf  to  a  maxwellian  distribution  at  the  initial  electron 
temperature 


CONCLUSIONS 

Collisional-radiative  model  coupled,  selfconsistently,  with  the  solution  of  Boltzmann  equation  has 
been  developed  for  hydrogen  atoms.  This  model  can  be  easily  extended  to  other  atomic  systems  and 
mixtures.  Results  have  shown  that,  in  particular  when  electron  temperatures  are  low  and  electron 
density  high,  eedf  differs  greatly  from  equilibrium.  Especially  in  these  cases,  time  dependent 
selfconsistent  calculations  are  necessary.  In  our  numerical  simulations,  we  have  tried  to  understand  the 
influence  of  the  variation  of  some  variables  (e.g.  gas  and  electron  temperatures)  over  level  populations 
and  eedf.  The  conditions  studied  are  a  rough  approximation  of  the  behaviour  of  a  plasma  produced  by 
laser  ablation  or  LIBS.  This  last  technique  pretends  to  give  quantitative  analysis  of  surface  composition 
and  plasma  temperature,  making  the  assumption  that  the  plasma  is  in  local  thermodynamic  equilibrium 
(LTE)  because  the  pressure  is  very  high  (1  atm).  In  this  work  we  have  shown  that  both  in  the  heating 
and  cooling  phase  the  system  is  far  to  be  in  equilibrium.  An  improvement  of  the  model  could  be  found 
in  a  self-coupling  with  a  fluid  dynamic  system,  which  should  give  more  informations  on  the  behaviours 
of  the  plasma  produced  by  the  interaction  with  a  laser  during  the  expansion  (see  for  example  [15,16]). 
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